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Abstract 



There is no known polynomial-time algorithm that can solve an NP problem. Evolutionary 
search has been shown to be a viable method of finding acceptable solutions within a reasonable 
time period. Recently quantum computers have surfaced as another alternative method. But these 
two methods use radically different philosophies for solving NP problems even though both search 
methods are non-deterministic. This paper uses instances of SAT, 3SAT and TSP to describe how 
these two methods differ in their approach to solving NP problems. 
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1 Introduction 



The mathematical study of the selection, grouping and permutation of a (usually) finite set of objects 
is called combinatorial optimization. Although real-world combinatorial problems have been known for 
centuries, it was not until the last 30 years — i.e., until computers because readily accessible — that finding 
a solution actually became practical. For example, the traveling salesman problem (TSP) requires one 
to specify the order in which a salesman should visit a fixed set of cities. The objective is to find that 
order which minimizes the total distance traveled. If N cites are to be visited, a brute force search must 
evaluate AH possible tours. Clearly, this is not practical for N much greater than eight or so. 

The real challenge of combinatorial optimization is to create algorithms and techniques that can solve 
realistically sized problems within a reasonable number of computational time JjJ. Most algorithms 
formulate a combinatorial problem as a search problem. Implicit is the idea that the solutions to 
combinatorial problems reside in an abstract solution space and two solutions are neighbors if they differ 
by a single mutation of a problem parameter. Associated with each solution is a real number that reflects 
fitness or quality of that solution. This space and the associated fitness values form a fitness landscape^. 
Any algorithm that "solves" an combinatorial problem is therefore a search algorithm that explores the 
fitness landscape. 

Unfortunately, many real-world combinatorial problems require such huge computational resources 
that brute force search methods are useless; they simply take too much time to find the optimal answer. 
This has led researchers to use search heuristics that yield an acceptable compromise: a possibly lower 
quality answer but with a minimal search effort. Evolutionary Computation (EC) techniques are at the 
forefront of this work and impressive results have been achieved. Nevertheless, these EC techniques still 
run on classical computers that use the Von Neumann model. But recently an entirely new approach 
has surfaced with potentially enormous consequences. This new approach is called quantum computing, 
which relies on the principles of quantum mechanics to evolve solutions. Existing programs, written in 
conventional high-level languages cannot run on these machines. In fact, no one even knows how to build 
one! Nevertheless, several system architectures have been proposed and algorithms are being developed, 
albeit in abstract form. 

It is interesting to compare how a quantum search, running on a quantum computer, differs from an 
evolutionary search, running on a classical computer. However, the whole point of this comparison is not 
to advocate one method over the other — its purpose is to highlight the radically different philosophical 
approaches. (Besides, because no one has ever built a quantum computer, there is no way any direct 
comparison can be made at this time. It is up to the reader to decide which approach holds the most 

1 In practice, fitness will be with respect to one or more attributes such as cost or power consumption; high fitness is 
associated with good values of the attribute. 
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promise.) If nothing else, the reader should come away with an appreciation for the total re-orientation 
in thinking that quantum search will require. 

The paper is organized as follows. Section ^| provides a broad overview of evolutionary algorithms and 
quantum computing. Because the focus is on NP-complete and NP-hard problems, a formal definition 
of these problem classes is also provided. Section ^| compares two quantum search methods against 
evolutionary algorithms for two well known NP-complete problems and one NP-hard problem. Finally, 
Section |] comments on the future of quantum computing. 

2 Background 

This section reviews evolutionary algorithms, quantum computing, and algorithms complexity. Each 
topic is wholely contained in a separate subsection so the reader may skip familiar material. 

2.1 Evolutionary Algorithms 

This section gives a brief introduction to evolutionary algorithms (EAs). For brevity, the emphasis is 
on those characteristics of specific interest to our work. More general information on EAs can be found 
elsewhere (e.g., see Q). 

Historically, evolutionary computations have a rich past, being independently developed by at least 
three independent research efforts, which ultimately produced three distinct paradigms: genetic algo- 
rithms, evolutionary programming and evolution strategies. All EAs share the same basic organization: 
iterations of competitive selection and random variation. More specifically, each generation (iteration of 
the EA) takes a population of individuals (potential solutions) and modifies the genetic material (prob- 
lem parameters) to produce new individuals via stochastic operations. Both the parents and offspring 
are evaluated but only the highest fit individuals (better solutions) survive over multiple generations. 
Although there are several varieties of EAs, they are all biologically inspired and generally follow the 
format depicted in Figure 111. 



initialize 
population 








evaluate 
fitness 








select 
survivors 








randomly vary 
individuals 









Figure 1: The canonical EA. Each individual in the population represents a unique solution to the 
optimization problem. The "select survivors" block chooses highly fit parents for reproduction while the 
"randomly vary individuals" block generates the offspring (i.e., new solutions for evaluation.) Selection 
may be deterministic or stochastic but generation of offspring is strictly stochastic. The process continues 
until either a fixed number of generations have been produced, or an acceptable solution has been found. 

The evaluation function for an individual returns a numeric value representing the quality of the 
solution described by that individual. This numeric value is often called the fitness of the individual 
while the evaluation function is called the fitness function. High fitness means the associated individual 
represents a good solution to the given problem. The selection process targets highly fit individuals for 
survival. 
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2.2 Quantum Computing 

Almost twenty years ago Richard Feynman observed that classical computers could not simulate certain 
quantum mechanical effects [[| . This observation spawned interest in the field of quantum computing — 
i.e., computational machines that perform calculations by emulating quantum mechanic effects. Al- 
though no practical quantum computer has yet been built (and the likelihood of building one in the near 
future is bleak), the interest in this held has not diminished because of the enormous computational 
potential such a machine can provide. Indeed, interest in this emerging held is growing by leaps and 
bounds [|. 

This section reviews some of the basic concepts of quantum mechanics that relate to quantum com- 
puting. No attempt has been made to be comprehensive and the reader is encouraged to study some 
of the excellent on-line tutorials j^, ^, |?j or journal reviews [||. (However, see the commentary on this 
topic in Section El.) 

Classical computer systems represent a single bit of information deterministically: the value is cither 
a logic or a logic 1. Quantum computer systems represent a single bit of information as a qubit, which 
is a unit vector in a complex Hilbert space C 2 . The ideas are commonly expressed using the bra/ket 
notation introduced by Dirac P]. The ket symbol is denoted by \x) and the corresponding bra is denoted 
by (x|. The ket describes a quantum state and the corresponding bra is its complex conjugate. 

In computer science domains the ket (bra) can be thought of as a column (row) vector. That is, the 
orthonormal basis {|0), |1)} can be expressed as {(0, 1) , (1, 0) }. Any complex linear combination of 
two kets is also a ket. The inner product of two vectors is denoted by (x\y). Note that since |0) and |1) 
are orthonormal, (0|1) = 0. \x)(y\ denotes the outer product of the vectors. 

Any practical quantum computer must manipulate a register of n qubits. If each qubit has an 
orthonormal basis {|0), |1)}, then a n qubit system has a basis expressed by the tensor product: C 2 <g> 
C 2 ® • • • C 2 . This gives 2™ basis vectors of the form 

|0)®|0)®--.|0) 

|0)®|0)®--.|1) 



|1)®|1)<8---|1> 

In general, \a) denotes the tensor product \a n ) ® |a„_i) • •• |ai) ® \ao), which means a quantum 
register has the value a = 2°a + 2 1 a i + • ■ ■ 2 n a n . 

A qubit need not exist in only one basis state. Indeed, a qubit can exist as a linear superposition of 
basis states co|0) + ci|l), where Co, c\ are complex numbers with |co| 2 + |ci| 2 = 1. More generally, the n 
qubit register can be prepared in a superposition of all possible classical states: 



(1) 



i=0 



where the normalization condition 



4 = i (2) 



must hold. The complex number Ci is called the amplitude associated with the state \i). 

The linear superposition of states is key to understanding how quantum computers operate. This 
linearity feature means that any operation on a superposition of states renders the superposition of 
that operation on each state individually JTo| ]. There is no analogue in classical computer system for 
this principle and, as will be shown below, it is an important ingredient of the power behind quantum 
computing. However, superposition also permits the following rather bizarre situation. Consider the 
state 

^(|00> + |11» 

This state cannot be expressed in terms of the individual qubit states. The proof is straightforward. Note 
that 

(<zi|0) + 6i|l>) ® (a 2 |0) + 6 2 |1>) = aia 2 |00) + ai6 2 |01) + M 2 |10) + &i6 2 |H) 

|00) + |11> 
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Clearly a\bi = 0, but this implies either a\ai or b\bi must equal zero, which is not possible. States 
that cannot be described by individual qubit states are called entangled. There is considerable debate 
concerning the actual role entanglement plays in search operations. This issue will be discussed in more 
depth in Section [| 

The state of a qubit register is determined by a measurement. In quantum systems this measurement 
process projects the system state onto one of the basis states. Referring to Eq. (|l|), the measurement 
returns a value of \i) with probability \ci\ 2 . Any subsequent measurement returns the state \i) with 
probability 1, which means the measurement process irreversibly alters the state of the system. Mea- 
surement also gives another perspective on entanglement: two qubits are entangled if and only if the 
measurement of one effects the state of the other. 

A quantum computer can perform the same function / as a classical computer if that function is 
a one-to-one mapping from the domain to the range. In other words, / must be a reversible function. 
Reversibility is usually mentioned in the context of performing computations without expending heat |y]] . 
Here, however, reversibility must hold or / will not be physically realizable on a quantum computer. 
Hogg ||Io|| illustrates the importance of reversibility with a simple example. Suppose /(si) = f(s2) — S3. 

Then for the superposition \s) = f- (|sx) + \ s 2)) linearity forces f{\s)) = — i — + l/( s 2))- But 

v( 2 ) v( 2 ) 
this equals y/\2)\ss), which violates the normalization condition given in Eq. (|J). 

Quantum systems evolve from state to state according to Schrodinger's equation E3|. Vector states 
can be expressed as a superposition of basis states each having an amplitude \ip%). This means evolution 
occurs by modification of the state amplitudes. Clearly, we would like to increase the amplitude of that 
state with the desired answer. Suppose we start in state |a) = ^] iph\o-k)- This system evolves over time 
under a linear operator U, i.e., \a') = U\a) = ^V^l^fe)- Hence, ip' = Uip and Eq. (Q) is satisfied iff U is 
unitary. To see this, consider the inner product (ip')'ijj' , which must equal one because state vectors are 
orthonormal. Substituting ip' — Uip yields 

{Uip)^{Uip) = ^(U^U)ip 

This inner product equals one if (U*U) — I. Hence, U must be unitary. 

It is convenient to adopt a simplified notation when describing unitary operations that are applied 
to individual qubits. Some common unitary operators are 

/: |0)-|0) 
|1>-|1) 



X: |0)-|1) 
|1>-|0) 

Z : |0> -> |0) 
|1>--|1) 

where / is an identity operator, X a negation operator, and Z a phase shift operator. Suppose we have 
a 3 qubit register and we want to negate the first qubit and leave the other qubits unaltered. This 
transformation is denoted by X ® I ® I. 

An extremely important transformation is the Walsh- Hadamard transformation defined as 

H: |0>-^(|0> + |1» 
|1>-^(|0>-|1» 

When applied to |0), a superposition state is created. When applied to n bits individually, a superposition 
of all 2 n states is created. Specifically, 

{H(g>H® •••<8>iJ)|000 ••■()) 
= -L((|l) + |Q))®(|l) + |0))®-.®(|l) + |0))) 
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y \x) 

On i-^i 1 ' 

Z x=0 

It is important to emphasize the role superposition plays in quantum computing. Let Uf be a 
unitary transformation corresponding to a classical function /, i.e., Ut : \x)\y) — > \x)\y © f(x)), where 
© represents bitwise exclusive-or. Notice that this transformation preserves the input — which must be 
done if / is not invertible — thereby making Uf unitary We can think of \y) as the hardware of the 
quantum computer. When this Uf operates on a superposition of states as in Eq. (|l|), the result is 

2™-l 2"-l 



Uf(y cMo)) = £ *u f (\i)\o)) 

8=0 8=0 

2™-l 

i=0 
2™-l 

= Eci|i)i/(»)) 



Notice that / is simultaneously applied to all basis vectors. Hence, a single application ofUf computes all 
2 n values o//(0), . . . , f(2 n — 1) at once jL4|. It is this quantum parallelism that is primarily responsible 
for the enormous interest in quantum computing. But something must be wrong. How can you extract 
an exponential amount of information out of a linear number of qubits? The answer lies with the 
amplitudes. If Cj = Cj Vi, j, then a measurement will produce any of the 2™ states with equal probability. 
Furthermore, once that measurement is taken, the system collapses into that measured state and all 
other information is lost. (You really can't get something for nothing.) Nevertheless, you can exploit 
this parallelism using the property of quantum interference. Interference allows the exponential number 
of computations performed in parallel to either cancel or enhance each other. Feynman [jl2| beautifily 
describes how light waves can constructively or destructively interfere to produce this effect. The goal of 
any quantum algorithm is to have a similar phenomena occur — i.e., interference increases the amplitudes 
of computational results we desire and decreases the amplitudes of the remaining results. It is a unitary 
operator that would alter these amplitudes. Examples of this approach are presented in Section |^. 



2.3 NP-Complete vs. NP-Hard Problems 

Many papers that discuss NP-complete and NP-hard problems (incorrectly) presume the reader fully 
understands the difference between these two labels. The distinction is important. For example, if one 
formulates a TSP problem as "does a tour exist of length < k" , then this is NP-complete. However, 
if the problem asks "what is the minimum length tour" , then this problem is NP-hard because it isn't 
in class NP. (They answer is not verifiable in polynomial time. The only way of answering 'yes' is 
to enumerate all possible tours.) In this paper a more formal approach is taken: all terms and three 
example problems are formally defined. This material is primarily taken from Jlq ], I begin with the 
following basic definitions: 

Definition 1 (decision problem) 

A problem in which the only answer is either YES or NO. 



Definition 2 (language) 

The set of all possible input strings to a decision problem that render a YES answer. 

The input strings are defined some fixed alphabet of symbols. For example, binary strings are defined 
over the alphabet {0, 1}. 

Definition 3 (polynomial-time algorithm) 

An algorithm that completes execution in a time which is a polynomial factor of the size of its input 
parameters. 



5 



Definition 4 (polynomially reducible) 

Let h\ and L 2 be two languages. L\ is polynomially reducible to L 2 (denoted by L\ oc L 2 ) if there exists 
some polynomial- algorithm that converts every input instance i% S L\ into another i 2 E Li- 
lt should be stressed that reducibility is asymmetric. In other words, if L\ oc L 2 , then this does not 
necessarily mean L 2 oc L\. Nevertheless, polynomial reducibility does have an important characteristic, 
which is given in the following theorem: 

Theorem 1 If Li oc L 2 , and there is a polynomial-time algorithm for L 2 , there there is a polynomial- 
time algorithm for L\. (See JZj|/ , page 343 for proof.) 

Definition 5 (nondeterministic algorithm) 

An algorithm that permits more than one possible move at some step during its execution. 

With these definitions it is now possible to define the algorithm classes P and NP. 
Definition 6 (class P) 

P is the class of languages (decision problems) L that, with input x, can in polynomial time return an 
answer YES if and only if x € L. 

Definition 7 (class NP) 

NP is the class of languages (decision problems) that can be checked for correctness in polynomial time. 

Notice that the above definition says nothing about the computational effort required to get that 
answer — it merely says to verify the correctness of an answer takes only polynomial time. Whether or 
not P=NP has yet to be determined. 

It is now possible to formally define NP-hard and NP-complete. It should be emphasized that the 
two type of problem classes are not interchangeable. 

Definition 8 (NP-hard) 

A problem V is NP-complete if every other problem in NP is polynomially reducible to V 
Definition 9 (NP-complete) 

A problem V is NP-complete if (1) V ENP, and (2) every other problem in NP is polynomially reducible 
to V 

NP-complete problems are decision problems. NP-hard problems ask for the optimal solution to 
an NP-complete problem. And, they have at least the same level of difficulty to solve as does the 
corresponding NP-complete problem. There are a very large number of problems that have been proven 
to be NP-complete 2 . Theorem 2 shows the most common way of proving a decision problem is NP- 
complete: 

Theorem 2 A problem X is an NP-complete problem if (1) X belongs to NP, and (2) Y is polynomially 
reducible to X, for some problem Y that is NP-complete. (See /pj[/ ; page 346 for proof). 

Theorem 2 is also used to prove if a problem X is NP-hard. Consider a problem X = (D,rj) that has 
an input domain D and some property r\. An algorithm which solves problem X uses an input instance 
I E D and verifies whether rj holds for input /. Suppose there exists another problem X 1 = (D',77') for 
which X oc X' . This means (1) I G D can be transformed to I' E D' in polynomial time, and (2) for 
any I E D, r] holds if and only if rf (based on I' E D') holds. X' will be NP-hard if X is NP-complete 
and X oc X'. 

I now describe two known NP-complete problems |l6j, which will be used in Section || Let S be 
a Boolean expression in Conjunctive Normal Form (CNF) — i.e., S is the logical and of two or more 
clauses where each clause is the logical or of Boolean variables or their complements. An example is 
S = (x + y + z) ■ (x + y) ■ (y + z). This Boolean expression is considered satisfied if an assignment of Os 
and Is to the Boolean variables makes S equal to a logic 1. 
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SATISFIABILITY PROBLEM (SAT): 

Instance: a Boolean expression in CNF 

Question: Does there exist an assignment of Os and Is to the variables such that the expression is 
satisfied? 

3SAT PROBLEM: 

Instance: a Boolean expression in CNF with each clause having exactly three variables 
Question: Does there exist an assignment of Os and Is to the variables such that the expression is 
satisfied? 

Finally, I describe a famous problem, which will also be used in Section 0. This problem is NP- 



TRAVELING SALESMAN PROBLEM (TSP): 

Instance: a finite set C = {ci,C2, . . . ,c m } of cities, and a distance d(ck,Cj) G Z + for each pair of 
cities Cfc, Cj S C. 

Question: What permutation [c n m, c w ( 2 ), . . . , c T ( m )] of C will minimize the tour length 



3 Search Approaches 

This section provides examples of how evolutionary and quantum search approaches have been applied to 
NP problems. No attempt has been made to survey the field; the objective is to present a few examples 
so the reader can appreciate the philosophical differences, which will be discussed in depth in Section ||. 

3.1 Evolutionary Search 

Although a number of papers have appeared discussing attacking SAT problems using EAs, I will focus 
on the recent work by Back, et al p7| . They used an evolution strategy to find solutions to instances of 
the 3 SAT problem. 

The search for a satisfiable solution is difficult because, as the authors point out, the fitness landscape 
is extremely flat — any genetic search reverts to a random search. Moreover, this type of landscape makes 
it difficult to define fitness in a meaningful way. The authors get around this situation by adapting a 
method suggested by a colleague This alternative method replaces each literal with x with (y — l) 2 
and x with (y + l) 2 . Furthermore, each disjunction A is replaced by an arithmetic + (sum) and each 
conjunction V is replaced by an arithmetic • (product). The resulting fitness function has a minimum of 
0, when the j/j's converge to 1 (true) or -1 (false). These changes convert 3SAT into a real-parameter 
optimization problem, which evolution strategies are ideally suited for. 

The evolution strategy randomly initialized the object parameters to values between -1.0 and 1.0. 
A (15,100)-ES with one standard deviation (<r) was used; a had an upper limit of 3.0. Later versions 
introduced various forms of recombination, which ultimately was shown to render the best performing 
version. 

3.2 Quantum Search 

Quantum search approaches differentiate between structured problems, where partial solutions can be 
extended to complete solutions, and unstructured problems. The unstructured approach can be used for 
finding solutions for NP-hard problems. 

2 A large database can be found at http://www.nada.kth.se/viggo/problemlist/compendium.html. 
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3.2.1 NP-Complete Problems 



Ohya and Masuda |19j developed a quantum search method that is frequently used for NP-complete 
problems. Their algorithm starts with the quantum system in the state 

I s ) = 4^ E ®?=ito>®?|0>®|0> 

for a SAT instance with variables Xx, . . . ,x n . The k qubits are garbage bits needed by reversible logic 
gates and the final qubit (initialized to |0)) indicates if the expression is satisfied. Then, using a unitary 
operator J7/, 



I*) = U f \s) 

i 1 

= ~m E ^ ®"=i ®* 1°) ® 1°) 



3?i ,. . . ,3:^—0 
1 



xi ,. . . ,x~ —0 



where /() is the Boolean expresssion. The last qubit is then measured by applying a projector P = 
I® |1) (1| to \t). If out of the 2" possible solutions there are r solutions that satisfy /, then the probability 
of measuring a solution is |P|t)| 2 = r/2™. For small r this probability is quite small. Hence, in practice 
quantum search algorithms try to exploit quantum interference to amplify the amplitude of the desirable 
solutions and attenuate all other amplitudes. 

Cerf, et al. |2Cf| provide a good description of exactly how this is done. Their approach relies on an 
"oracle" function f[x) that equals one for the optimal input x (and zero elsewhere). The goal is for the 
quantum system to evolve from an initial state \s) to the target state \t) in minimum time. Note that 
f(x) = 1 only at x = t. More precisely, the goal is to increase the amplitude of \t) to a point where a 
measurement will render \t) with the highest probability. 

Assume an arbitrary unitary operator U has been found that connects \s) to \t) — i.e., (f |f|s) ^ 0. The 
probability \t) is actually found is |(i|[/|s)| 2 , which means the experiment must be repeated |(i|t/|s}| -2 
times on average to guarantee success. However, it is possible to reduce this to the order of | (t\U\s) | _1 — 
which can be a considerable savings — with an appropriate quantum search algorithm. 

The algorithm begins in a superposition of states and any measurement is postponed until the end. 
Cerf, et al. defined a specific unitary operator 

Q = -Ue™ p >tfe™ p * 

where P s = \s)(s\ and P t — \t)(t\ are projection operators on \s) and \t). These exponential operators 
simple flip the phase on a state. For example, the phase of state \x) is flipped by e i7rPs iff x = s. Since 
the objective is search for state \t), the oracle is used to implement its exponential operator. That is, 
e l7rPt \x) = (— iy^\x). Then, by repeatedly applying Q, the amplitude of \t) is amplified, beginning at 
U\s). This amplitude amplification is achieved by the repeated application of Q which, in effect, rotates 
the starting state \s) into the target state \t). In other words, the beginning state U\s) is rotated to the 
target state |t) by repeated applications of Q, followed by a measurement. Recall U was an arbitrary 
unitary operator; by using structure information it may be possible to find a better U 1 so that U'\s) has 
larger amplitudes in states which are more probable to be solutions. Cerf, et al. p3] describe a method 
that constructs such a U'. 

Grover's quantum search algorithm searches a random database of N items in 0(t/~(N)) steps [pl| . 
This means unstructured NP-complete problems can be solved by forming a database of all possible 
candidate solutions, and then use Grover's algorithm to find the solution. Although this is a considerable 
speedup over classical machines, it may not be all that impressive. For instance, if one has to find an 
assignment of one of k values to n total variables, a classical algorithm would take 0(k n ) steps while 
quantum algorithms would still take 0(k n ^ 2 ) steps. Nevertheless, the algorithm does find a use with 
both NP-complete and NP-hard problems. 
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3.3 NP-Hard Problems 



A beautiful example of how non-traditional architectures can solve NP-hard problems is the scheme 
presented by Cerny p2| to solve an instance of TSP. Figure || shows an interference experiment setup. 




P2 = I(|) 2 I 2 



Figure 2: An interference experiment showing how amplitudes combine in both a constructive and 
destructive manner. The electron gun fires electrons that go through slits in a wall. A movable detector 
determines where the electrons impact the backstop. The wavelike behavior of the electrons produces 
interference so that the total distribution P12 7^ Pi + p2- This figure was adopted from Jl^ |. 

Cerny proposed a quantum computer similar to that of the interference experiment. This computer 
has (n-1) walls representing cites 2,3,. . . , n. Furthermore, each wall has (n-1) slits. A beam of quanta 
(e.g., electrons) sent through this array has possible trajectories. The wavelike behavior of 

electrons means a superposition of all possible trajectories is rendered in 0(n) time. A sample trajectory 
in this quantum computer is shown in Figure [|. 
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Figure 3: A sample TSP trajectory in the Cerny quantum computer. This tour is S, 3, 2, 5, 4, D. Note 
that not all tours are "legal" . For instance, S, 3, 3, 5, 4, D is also a trajectory but it is illegal because city 
3 is visited twice and city 2 is never visited. 

These trajectories identify tours but they do not indicate the length of those tours. Since this 
machine is hypothetical, an internal degree of freedom can be added — even if nature doesn't provides it. 
Specifically, the internal state of a particle is 

|fc;c 2 ,c 3 , . . . ,c„;p) 

where k £ {0, 1, 2, . . . , NL}, Ci S {0,1}, and p S {0,1}. The quantum number k measures the tour 
length; Ci = if city i is not visited and 1 otherwise; and the quantum number p is used to control the 
dynamics of the search ^. 

3 See the appendix in |^2[] for an explanation of how p is used; its purpose is not needed for the brief overview given in 
this paper. 
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To illustrate the dynamics, let (i, m) — > (i + l,n) denote a trajectory between slots on neighboring 
layers i and i + 1 indicating the tour moves from city m to city n. If the particle moves through slot 
(i, n), then 

c n > c n — 1 . 

Furthermore, assume the particle moving between layers encounters a field that increases the quantum 
number k by a factor d mn if the trajectory moves from (i, m) to (i + l,n), where <i mn is the distance 
between the two cites. Then, with an initial state 

|0;0,0,...,0;0) 

, after passing through the machine the particles are in a state 

1^' C2 > c 3i ■ ■ ■ > c ™j ^trajectory- 

trajectories 

Note that the legal tours will have all Ci = 1 and the quantum number k is the tour length. A filter 
installed at point D purges all kets with at least one a = 0. This leaves 

|fc;c 2 ,c 3 ,...,c„;p) TS route . (3) 

TS routes 

A Stern- Gerlach-like device Q could be used to construct such a filter, which would split the above 
superposition into NL streams according to the k value. A set of particle detectors would then indicate 
the tour length — i.e., a detector measuring M would fire if there exists a TS tour with a length equal to 
M. 

It is important to emphasize this does not mean an instance of TSP can be solved in polynomial 
time. In principle, Eq. (^) could represent a superposition of 0{n\) states. Hence, even if Grover's 
algorithm is used, it would take 0{y n\) steps to find the minimal length tour. 



4 Discussion 

No one has yet built a quantum computer capable of searching for solutions to even moderate size NP 
problems. But, despite our inability to make head-to-head comparisons of evolutionary and quantum 
searches, it is possible to highlight their primary philosophical difference: 

evolutionary search 

The algorithm uses stochastic operations to explore a fitness landscape comprised of all possible 
solutions. 

quantum search 

The algorithm forms a superposition of all possible problem states and then applies a unitary 
operator to compute a superposition of all possible solutions. This unitary operator also alters the 
amplitudes of each state exploiting interference to maximize the amplitudes of the desired states. 
A final measurement extracts the solution with a probability equal to the amplitude squared. 

EAs must tradeoff exploration against exploitation. In other words, the EA must carefully decide 
which regions of the fitness landscape to abandon, because the solutions are found there are poor, without 
putting much emphasis on regions with good solutions because that would limit the search. The focus of 
EA research with respect to NP problems is in two areas: (1) identification of appropriate representations 
of the problem parameters, which ultimately defines the fitness landscape, and (2) creation of effective 
stochastic reproduction operators that control movement over the fitness landscape. 

Quantum search algorithm exploit superposition to produce massive parallelism. One rather con- 
tentious debate in this field is the role entanglement plays. On one side of the fence are those who 
feel entanglement is essential for speedup |23fl , while on the other side are those who feel it is com- 
pletely unnecessary [^4| |25|. The latter group believes superposition and interference are sufficient to 
produce speedup. This issue could be resolved if a truly entangled system were available for study. 
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Unfortunately, recent room-temperature liquid-state NMR experiments have failed to produce any en- 
tangled states. Still, some researchers feel increasing the number of qubits (currently only around 3) will 
eventually make entanglement appear p6[ . 

One of the main difficulties in running a quantum computer is they must remain completely isolated 
from their environment or the state evolution will cease. Furthermore, there is no way of observing 
what's going on unless a measurement is taken. But taking a measurement process changes the system 
by causing it to collapse into one of the basis states. Some methods of dealing with this have been 
proposed ||, but it still remains a thorny issue. Consequently, we can expect running a quantum search 
will be much more fragile than the running of an evolutionary search on a classical computer. 

One final note on unstructured NP-complete problems. The 0(-y/{N)) time for Grover's search algo- 
rithm has been proven to be optimal [^J. This has a rather disappointing consequence: if 0(-\/{N)) time 
is optimal, this may mean quantum computers can't solve NP problems with an exponential speedup. 
Preskill p8| suggests the real application area may lay outside NP. Quantum system simulation is one 
example, which was also previously suggested by Feynman 0. 



5 Final Comments 

I will conclude this paper with some personal observations. I do believe quantum computing will change 
the way computer engineers and scientists think about computing systems. To date, quantum computing 
has been the domain of primarily physicists. It is about time that computer engineers and scientists 
enter this arena and begin to drive its direction. 

Many computer professionals entering this field are quickly put off by the notational and conceptual 
barriers. Tutorials are available (e.g., see |Q |2^, ^o|, ||), but many readers will quickly find them 
incomprehensible — they are written by physicists for physicists. (Out of this lot, however, I believe [^) is 
the best.) The sad truth is a computer professional who lacks a firm foundation — i.e., formal training — 
in quantum theory will most likely not be able to contribute to the quantum computing field. As an 
absolute minimum I would recommend an upper division undergraduate course in quantum mechanics. 
This should be sufficient background for one to begin reading the literature from the field. 

My other observation concerns the practicality of the currently proposed quantum computer architec- 
tures. Many proposed systems (including the NMR approaches) contain a vast network of interconnected 
quantum gates (such as AND gates) which implement some function f(x) I believe this is entirely 
too low of a design level, which is unlikely to lead to massive improvements in computation power — 
certainly no where near orders of magnitude improvement. Although, in principle, all computer systems 
are just interconnected primitive logic gates, engineers typically do not visualize them in this way. For 
example, few designers think of a processor as a network of primitive logic gates implementing Boolean 
expressions. For the most part architectural design work is rarely performed lower than the register- 
transfer level. Indeed, the hardware description languages in use today, such as VHDL and Verilog, are 
most frequently used at the register-transfer level. 

Thinking of quantum computers in terms of interconnected logic gates also tends to limit their ability 
to perform general purpose computations — especially those computations that are inherently parallel. 
For instance, can a quantum computer perform an evolutionary search? 

I am convinced that a radical increase in computing power will only come once the Von Neumann 
paradigm has been dispensed with. Architectures such as those proposed by Cerny |]22f are an example 
of the imagination that will be required. 
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